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Abstract 


The effects of particle drift as well as diffusion and convection are 
employed to describe the propagation of solar charged particles in the inner 
solar system. It Is found that protons and electrons gradient drift in 
opposite polar directions in the interplanetary magnetic field, and al- 
though the charge separation thereby produced is very small, this drift 
determines the energy loss of these particles. In addition, it is shown 
that electric field drift causes protons and electrons, alike, to co- 
rotate but slower than the magnetic field structure. As a result, solar 
particles injected onto a given magnetic field line sweep past the earth 
later than the field line. The resulting azimuthal intensity profile is 
determined by the distribution of sun-earth transit times for these par- 
ticles as they diffuse along the magnetic field lines. These considerations 
are applied to low energy protons and electrons and predictions for their 
energy loss and intensity profiles are made. A comparison of these pre- 
dictions with the data from the Imp 1 satellite indicates that the observed 
recurrent protons might be understood in terms of this analysis if the dif- 
fusion is sufficiently slow to give 1 Mev protons a 4 day transit time. 

The observed electrons, however, cannot have had such a solar origin and 
are evidently galactic. 
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I. Introduction 

In order to correctly understand the properties of galactic cosmic 
rays and particles of solar origin it is necessary to understand the 
effects produced on these particles by passing through the interplanetary 
field and solar plasma. This problem has, of course, been extensively 
studied for years and has been reviewed in detail by Parker (1963)* 

Ouring recent years much has been learned about the solar wind and 
the interplanetary field, and this has now made possible detailed studies 
of the propagation of charged particles through the solar system. In 
particular during 1963*1964/ the Imp I satellite provided extremely valu- 
able data by carrying detectors which simultaneously made measurements on 
the interplanetary magnetic field (Wilcox and Hess, 1965 )/ the solar plasma 
(Brld^gvet al, 2965 } low energy protons (Fan, et al, 19 66 ) and electrons 
(Cline, et al, 1964) and galactic cosmic rays (Balasubrahmanyan, et al, 1963 )* 
The data from Imp I which is of interest to us in the present work is shown 
in Figure 1. At that time the i nterplanetary magnetic field formed four 
sectors of alternating polarity whose boundaries are indicated in the 
figure by vertical lines. (These lines appear cross hatched when the exact 
time of passage of the sector boundary is uncertain due to the satellites 
being near perigee.) As one can see, low energy protons, which evidently 
do not originate from observable solar flares (Fan, et al, I966), were 
observed to recur in a particular positive magnetic sector while electrons 
of comparable energy seemed to recur predominantly in an adjacent negative 
magnetic sector. Also shown in the figure is the counting rate for galactic 
cosmic rays above 30 Mev/nucleon and the bulk velocity of the solar wind. 
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The observed association of recurrent particle beams with magnetic 
sectors of fixed polarity lead us to examine the possibility that these 
Mev particles are released from the sun preferentially along a magnetic 
sector boundary and subsequently undergo a drift from that boundary. 
Sections II and III are therefore devoted to the description of particle 
drift in the interplanetary magnetic field and its relationship to the 
transport of particles between the sun and earth. Drift has ordinarily 
been discounted as being small for solar flare particles (Parker, 1963) 
(Axford, 1965); however we find that drift is of great significance in 
determining the energy loss as well as perhaps, the intensity-time profiles 
for events of the kind described above. 

In the last section we investigate the possibility that the observed 
electron and cosmic ray time variations are simply the result of modula- 
tion by a variable solar wind. This has been done by evaluating the linear 
correlation coefficient between the daily particle intensity records and 
the plasma velocity records, and we find that cosmic rays, in particular, 
tend to follow the daily solar wind variations. 
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II. Drift in the Interplanetary Magnetic Field 


We shall refer to a spherical coordinate system (r,D,$) centered at 
the sun and at rest with respect to the earth and begin by considering the 
motion in the polar direction of charged particles ejected from a particular 
location on the solar surface. Assuming that the interplanetary 

magnetic field is reasonably regular in the inner solar system, the motion 
of protons and electrons of **10 Mev may be treated in the guiding center 
approximation. The most general drift velocity, to first order in the ratio 
of the Larmor radius to the field scale size, is given by Northrop (1935) 
as 


| ^ + ? *i»i +? <».df 

where M is the particle's orbital magnetic moment, tx is the unit vector 

cExB 

along the magnetic field and = -jjf- . Approximating that the fields 
are static, the drift velocity becomes 




where the A • V term, which is small, has been neglected. 

The interplanetary magnetic field is assumed to have the underlying 
Archimedes spiral pattern (Parker 19$) 
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where B is the field at the sun, r is the solar radius, SL is the solar 
o o 

angular velocity and is the solar wind speed. Remembering that the 

solar wind is radial and that £ * - 7 V x B, equation 2 becomes 
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In obtaining equation 4, it has been approximated that 8 is independent of 
( p over the region of interest and that the phenomena at hand occur near 
the ecliptic plane ( 0 ^ £ ) . The polar drift velocity becomes 
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Now the magnetic moment is given by and the guiding center velocity 

. B ^2 2 -1/2 

(which is approximately v^ ) has the components "r-18 l V » (1 
B a2 p “1/2 ~ 
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Then by use of equation 3> one 


obtains 
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We are interested in this drift velocity for r AU, corresponding to 
jlr£V o< and in this case the second term may be neglected. Also we expect 
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that magnetic irregularities in the field maintain a random pitch angle 
distribution for the particles, implying that on the average, Vj_ 2 * 2V ^ 
We therefore approximate 
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where is the magnitude and polarity of the magnetic field at 1 AU. This 
implies that protons and electrons of the same kinetic energy drift in the 
polar direction with equal speeds but in the opposite direction. If the 

magnetic field is outward (B^X)) equation 7 implies that protons drift 
south and electrons drift north, whereas the drift directions are reversed 
if the magnetic field is inward toward the sun. Therefore if protons and 
electrons are ejected from, say, the northern solar hemisphere, protons 
drift toward the ecliptic (and hence the earth) if the magnetic field is 
outward and electrons drift toward the ecliptic if the magnetic field is 
inward. It must be noted at this point that any drift in the polar direc- 
tion by these particles implies a change of energy by an amount 
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where B^ is in gauss, R^ in cm, and XL in sec . Defining as the 

particle energy when at 0 , F (0) * c - c 0 \ where 0*0-0 
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and c * 10“ s e il-B e R e 2 (1 + — — ~) . From equation 7 the angular drift 
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where V^(r > 0 / ) is the mean speed with which the particles move outward at 
r and0. Furthermore it follows from equations 6 and 8 that 
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This energy loss function has the limits, 
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and tells us that particles near the sun lose energy more slowly than has 
previously been thought (Parker, 1965 )* 

We now consider the drift of particles in the azimuthal direction. 
From equation 4, 
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implying that electric drift causes the particles to corotate with the field 
structure near the sun, but slower than the field structure with increasing 
r. fislatl ve to the field structure , the azimuthal drift velocity becomes 


do) 
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III. Bulk Motion in the Interplanetary Field 


A unified description of the simultaneous drift, convection and diffu- 
sion of particles propagating through the interplanetary medium can be ob- 
tained by constructing a Fokker-Planck transport equation for the problem. 
Letting N(r, £,t) represent the number of particles/volume with energy £ 
at the time t, and letting j represent the particle current density, the 
Fokker-Planck equation is 


3 ...r 


Tt* V ’ 1 *7e <N£) • ° ■ 


D iff us ion ; convection and drift all contribute to j so that 


~r j 

^ " -'diffusion * ^convection + ^drlft 


Since the scale of the interplanetary field structure seems to be large 
compared with the Larmor radii we consider here, the particle motion is 
well approximated as a guiding along the force Uines# But small scale 
irregularities in the magnetic field can be expected to produce non-ad iabatic 
changes in the pitch angle (Roelof, 1 9 ^ 6 ) and thus a particle tends to find 
itself with a mirroring pitch angle after some characteristic distance 
(the mean free path). This gives rise to a random walk along the lines 
of force and leads to the usual diffusion current, 
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The radial expansion of the solar wind convects the magnetic irregularities 
radially outward, but since the impulse given to the particles by the ir- 
regularities during the mirroring is along the line of force, the result- 
ing particle convection velocity is the component of the solar wind velocity 
along the field direction, or 
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where b is a unit vector along the local field direction. Finally there 
is a macroscopic current due to the particle drift normal to the field 
direction. Since we are here interested in the mean drift velocity of 
all particles, rather than guiding centers, in a unit volume, we must use 
the current density 
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rather than equation 1, {see Spitzer ( 1962 ))* As before, E « - x B. 
The Fokker-Planck equation then becomes 
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where V^.^. and £ are given in equations 9a and 13c. It seems worth 
pointing out that equation 13b and 13c imply that in the absence of dif- 
fusion, particles stream radially outward from the sun at large r, but 
that this is due to electric field drift and not convection. 

Rather than attempting to solve equation 14, we shall adopt an approxi- , 
mate treatment of the combined diffusion and drift which essentially amounts 
to first solving the diffusion problem without drift and then calculating 
the drift of the particles as though they were propagating outward at 

their diffusion velocities. Near the sun (SLr «V ) we have that 

o 

A • 

j e N V b, |V, ,, I « V and c is small. Furthermore, approximating 
J conv o * ' drift o 

that the diffusion is along the r direction and ignoring the divergence 
of space with increasing r, equation 14 lead simply to 
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Following Parker (1965) the solution of equation 14a corresponding to the 
boundary conditions that 


N(<?,t ,£)«0 and N(r, 0,£) 
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By differentiation with respect to time It follows that (for r .X << Kt) t m , 
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the time at which most particles with K arrive at r, is 


3K + (5K 2 + V 2 r®) ,/2 

t n « v*° (,5a) 


and therefore a diffusion velocity may be defined as 


¥ d dt 
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for these particles. Particles arriving at r in the time t have diffused 

t 

more rapidly by the factor, and may therefore be characterized by an 
average diffusion velocity. 
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We have therefore replaced the statistical diffusion of the particles by 
a steady outward propagation at the diffusion velocity in this approx ima- 
tion. It is now possible to estimate the drift of the diffusing particles, 
and according to equation 12, the azimuthal drift becomes 
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Consequently the above change of independent variable from t to ^ in 

equation 13 gives the azimuthal intensity distribution, N(R , 0 , £), to 

6 



( 21 ) 


- 1 6 - 


o 


c©' 



•39 6„, for 5 K«V c R e 

The diffusion coefficient may be expressed as 

( 10 10 X for relativistic electrons 
K = J 

^ .Vj x 10 s X£ q ^ 2 for non-relat ivistic protons 

and choosing X = 2 x 10 11 cm, = 300 km/sec, the final energy and 

the drift, Q'iQ may be evaluated. These functions are shown in Figure 2. 
The flat relative energy loss for the electrons and low energy protons is 
a reflection of the fact that these particles have transit times which are 
energy independent, and the protons, having longer transit times, lose a 
larger fraction of their energy. The actual amount of drift is seen to be 
small for all particles in the energy range of interest. Therefore although 
this drift could in principle produce a "charge separation" of protons and 
electrons, the angular separation is very small under the prevailing con- 


d i t ions. 
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We next evaluate the azimuthal drift of the particles* According to 
equations 12 and 19 * the location of the particles relative to the field 
structure is found from 



for 3*> V R 
o e 


( 22 ) 


Integrating, one obtains 
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which has the value 

(j) « *5^ radians for protons 

* *02 radians for electrons 

at £ = 1 Mev. That is, if 1 Mev protons and electrons are ejected from 

the sun onto a given line of force, one should expect that the protons will 
have drifted some 20° to the east of the line of force by the time they reach 
1 AU whereas the electron will stay essentially tied to the line of force. 

The form of the azimuthal particle intensity distribution follows from 
equations 15 and 17- Letting N(A',K) represent the intensity of particles 
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with diffusion coefficient K arriving at the azimuth 0 ', we have 


%SK) = N(t,K) ) 


( 24 ) 


which is evaluated for 1 Mev and 30 Mev protons in Figure 3 , with n(£) ■ £ - 4. 
The peak of the distribution for 1 Mev protons is understandably at a larger 
angle than for the 30 Mev protons since they drift azimuthal ly during a 
longer transit time. (It should be recalled at this point that, contrary 
to the polar drift velocity, the azimuthal drift velocity, equation l I, 
does not depend upon particle energy. The increased drift for the 1 Mev 
protons Is thus a reflection of the longer transit time only.) 


V. Azimuthal Variation of Energy Spectrum 


We now restrict ourselves to protons and examine the effects of pro- 
pagation on the energy spectrum. The variation of the energy spectrum with 
azimuth can be studied by obtaining the time variation of the spectral slope 

at some fixed energy. Referring to equations !$Ta and J? 3 , we find that 
t 00 

■ * " is very nearly energy independent. Thus, confining ourselves to 

^ m 

times and energies for which Kt » R X, equation 2 k becomes 
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where K - .87 x 10 £ and t - • Transforming to logarithms 
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we obtain 
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This time variation is displayed in Figure k at three energies. (In Figure 4 
and elsewhere we use three variables which must be distinguished: t is a 
transit time for particles passing between the sun and R ■ l AU, p' is an 
azimuthal angle measured from the magnetic sector boundary along which the 
particles are injected, and T ** ^ P 1 * s the time elapsed since the sector 
boundary swept past the earth.) In Figure 5 # curve a, the energy spectrum 
at f ■ 3 days is plotted as an example of a propagation-modulated spectrum. 
The following conclusions may be drawn from these curves. 
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(0 Only d ifferences in slope at different energies or times are of signi- 
ficance since the spectral index of the injection spectrum appears addi- 
tively in equation 25 . (2) The spectrum has the characteristic behavior 

of falling rapidly to a maximum steepness after approximately 2 days and 
then flattening very gradually with increasing time. This time dependence 
reflects the fact that high energy particles, having the shortest transit 
times, are found closest to the sector boundary while lower energy 

particles predominate more and more at greater distances from the sector 
boundary. At the solar wind transit time of 5*7 days 2 days) the 

maximum intensity of low energy particles arrives (with the solar wind) 
and the energy spectrum is therefore steepest. For larger values of T , 
the plasma carrying the low energy particles is beyond 1 AU and one is 
observing particles which have diffused inward, in the rest frame of the 
plasma. Therefore as increases one detects particles which have retro- 
gressed progressively farther and farther from the plasma and these are 
the higher energy particles. The energy spectrum therefore flattens again* 

(3) The time variation described above becomes less pronounced at higher 
dV 

energies because (« “ h° w ra P*^Y the transit time varies with 

energy, decreases with increasing energy. Therefore at high energies, where 
all particles have nearer the same transit times, the particles tend to main- 
tain their initial energy spectrum as they propagate outward. 

In this discussion, as in equation 2k, the variation of the spectrum 
due to deceleration has been ignored. Without solving for the energy spec- 
trum resulting from the coupled diffusion and energy loss in the present 
calculation, we simply indicate the effect of deceleration alone on the 
spectrum. If we suppose that the energy spectrum at the sun is n(€) « 


-4 
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and all particles arrive simultaneously at 1 AU with energies ■ g(£^ 
given by equation 21 then their energy spectrum becomes 

rn'lfy - "tfo* "iT 

This spectrum is plotted in Figure 5, curve b, and it may be seen that 
the deceleration is sufficiently small so that the spectrum Is modulated 
very little. We conclude therefore that the energy spectrum observed at 
the earth should rather resemble the energy spectrum at injection except 
at early times (fVO) and very low energies (€f il Mev). 
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VI. Application to Imp I Data 

We now turn to the data obtained from the Imp I satellite collected in 
Figure 1. The integral intensity profile for Mev protons is shown at the 
bottom and these recurrent shapes are to be compared with the predicted 
profile. Figure 3. The general features of a steeper ascent than descent 
and a peak intensity about two days after the boundary seem to be reflected 
in the data. Furthermore, Fan et al, (1966) report that the energy spectrum 
is slightly flatter (E~ 4 ) during the first days of the sector than in the 
last days (E* 5 ), and this is consistent with the predicted slope variation. 
Figure 4. Therefore it seems that if the sun released protons only along 
the sector boundary observed on 4 December and was inactive elsewhere, then 
the subsequent diffusion and drift of these particles would give rise to 
intensity profiles and energy spectra similar to those observed. 

The observed intensity profile for Mev electrons indicates that maxima 
occur three days or more beyond the 12 December sector boundary. However, 
our analysis predicts that solar electrons released from the sun along 
that sector boundary should be observed predominantly at the boundary, and 
therefore we conclude that the observed electrons have not had such an 
origin. 

Let us now consider the alternative idea (Brunstein, 1965) that the 
electrons are galactic. In Figure l the variation of the solar wind velo- 
city across the field sectors is displayed. Let us assume that this vari- 
ation persists for a distance AR beyond the earth. Then if the inward dif- 
fusion of cosmic rays is largely confined to the spiral field direction so 
that the cosmic ray density along a given line of force depends only upon 
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the local wind velocity, the modulation, given by Parker (1965)* is 


n E 


n o exp 


.2M (1 + l ( JSd R)2) 

U 3 * V ' ' 


(26) 


where n^ is the flux at the earth, n Q is the flux beyond and V is the 
variable wind speed. Choosing the values A * 2 x 10 iA cm, m l AU, 
c and SLA R *» 1»00 km/sec, one finds that the percent change in the 
galactic intensity corresponding to a wind velocity changedV should be 


dn F 

— - «*3x 10^ dV (27) 

n E 

However the observed percentage changes indicated in Figure 1 
are much larger than one would expect from equation 27* which predicts that 
a solar wind variation of 100 km/sec should produce only a 3^ change in the 
electron intensity. In order to quantitatively examine the extent to which 
the daily variations obey equation 2 6, we have studied the daily correlation 
between the electron intensity and the plasma velocity. In Figure 6 the 
correlation coefficient for various detectors is plotted as a function of 
delay time between the particle intensity and plasma velocity records, and 
it may be seen that the electron intensity shows no significant correlation 
with the plasma velocity for any value of the delay. This suggests that the 
modulation described by equation 2 6 cannot alone describe the propagation of 
the observed electrons. 

If equation 26 is applied to 30 Mev protons, having cf .3, one 


obtains 
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drip 

— - - 10" 8 dV (28) 

n E 

assuming the same parameter values as for equation 27* Thus a solar wind 
variation of 100 km/sec is predicted to produce a 10$ variation in the 
cosmic ray intensity, and the data does indeed indicate such variations. 

The daily correlation between the cosmic ray intensity and the plasma 
velocity is indicated in Figure 6 for each of four Imp I detectors (des- 
cribed by Balasubrahmanyan et al, 19^5) and the Deep River neutron monitor. 

In each case the data reveals a statistically significant correlation which 
is strongest when the cosmic ray data is compared with the plasma data of 
the previous day. However there is weaker, yet statistically significant, 
correlation for time delays ranging f rom ** +1 to -5 days between the 
cosmic ray and plasma velocity records. Without attempting to provide a 
detailed interpretat ion of these time delays in this paper, we only comment 
that negative delays are to be expected since the electric field drift 
(section 111) causes the particles to lag slightly behind the field lines 
about which they are guiding. Such a field line, characterized by a certain 
plasma velocity, therefore crosses the earth before the cosmic rays asso- 
ciated with that field line, in agreement with the negative delay in Figure 
6. The only point we wish to make at this time, however, is that the cosmic 
rays follow the transient behavior of the plasma velocity while the electrons 
do not. Since neither the model based on solar origin described earlier, 
nor the modulation model, equation 26, explain these electrons, it would 
seem that perhaps processes are operating on electrons which have little 
effect on the heavier cosmic rays. 
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Summary 

The interplanetary electromagnetic field, in a reference frame at rest 
with respect to the earth, has been represented by an Archimedes spiral 
magnetic field and the electric field induced by the solar wind. It is 
found that charged particles ejected from the sun drift in these fields 
in such a way that they tend to co-rotate with the magnetic field while 
being decelerated by the electric field. At the same time the particles 
diffuse along the magnetic field lines and arrive at 1 AU with a dispersion 
in transit times which in turn implies a dispersion in the drift. Making 
the simplifying approximation that the diffusion is one dimensional with 
a constant mean free path equal to ^.01 AU, proton intensity profiles 
and energy spectra have been calculated, and it is found that these pre- 
dictions are in reasonable agreement with the recurrent proton data from 
the Imp I satellite. Also, the cosmic ray data from Imp l was found to 
correlate well with the daily variations in the solar wind velocity. This 
implies that a '*stream" of low energy cosmic rays is essentially always 
in the presence of one and the same plasma velocity and hence, one ami the 
same bundle of field lines. 
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Figure Captions 


Figure 1. Imp I data , November, 1963-January, 1964* The vertical lines 
indicate magnetic sector boundaries and the crosshatching represents the 
interval of uncertainty in these boundaries* The first sector was charac- 
terized by an outward directed field and the succeeding sectors alternate 
in polarity* 

Figure 2. Energy loss and polar drift of solar charged particles* This 
energy loss is caused by the gradient drift of particles against the in- 
duced electric field of the solar wind. 

Figure 3- Predicted intens ity-t ime profiles for solar protons* The in- 
tensity is plotted against time measured from sector boundary passage* 

The corresponding angular position is also indicated* 

Figure 4* Predicted time variation of the spectral index for solar protons* 
The spectral index at each energy is plotted against time (and angular 
position) as in Figure 3» 

Figure 5* Representative differential energy spectrum for solar protons* 
Curve a represents the energy spectrum at T *= 3 days resulting from pro- 
pagation without energy loss for an injection spectrum E^”*. Curve b 
represents the energy spectrum resulting from energy loss only, for the 
injection spectrum, E q 


-4 


Figure 6. The time correlation between observed particle intensity and 
solar wind velocity* The linear correlation coefficient is plotted against 
delay time between particle intensity and plasma velocity records. Posi- 
tive delay indicates that the particle intensity was correlated with the 
plasma velocity measured at an earlier time. Iff , C (plastic scin- 

tillator) and omnidirectional intensity and electrons (3*12 ttev) represent 
measurements made with the Imp I satellite. 
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